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ABSTRACT 

The aim of this paper is to cast light on the formation and evolution of elliptical 
galaxies by means of N-body/hydro-dynamical simulations that include star forma- 
tion, feed-back and chemical evolution. Particular attention is paid to the case of 
dwarf spheroidals of the Local Group which, thanks to their proximity and modern 
ground-based and space instrumentation, can be resolved into single stars so that 
independent determinations of their age and star formation history can be derived. 
Indeed, the analysis of the color-magnitude diagram of their stellar content allows 
us to infer the past history of star formation and chemical enrichment thus setting 
important constraints on galactic models. Dwarf galaxies are known to exhibit com- 
plicated histories of star formation ranging from a single very old episode to a series 
of bursts over most of the Hubble time. By understanding the physical process driving 
star formation in these objects, we might be able to infer the mechanism governing 
star formation in more massive elliptical galaxies. Given these premises, we start from 
virialized haloes of dark matter, and follow the infall of gas into the potential wells and 
the formation of stars. We find that in objects of the same total mass, different star 
formation histories are possible, if the collapse phase started at different initial densi- 
ties. We predict the final structure of dwarf spheroidal galaxies, their kinematics, their 
large scale distribution of gas and stars, and their detailed histories of the star for- 
mation and metal enrichment. Using a population synthesis technique, star formation 
and metal enrichment rates are then adopted to generate the present color-magnitude 
diagrams of the stellar populations hosted by dwarf spheroidal galaxies. The simula- 
tions are made assuming the red-shift of galaxy formation Zf or = 5 and varying the 
cosmological parameters Hq and qg. The resulting color-magnitude diagrams are then 
compared with the observational ones for some dwarf spheroidals of the Local Group. 

Key words: Galaxy Formation - N-body Simulations - Stellar Populations - CM- 
Diagrams 



1 INTRODUCTION 



properties that are more compatible with very old stellar 
activity? 



The present-day challenge with elliptical galaxies (EG's) is 
to unravel their formation and evolutionary history. In a 
simplified picture of the issue, the problem can be cast as 
follows. Do EG's form by hierarchical merging of pre-existing 
sub-structures (maybe disc galaxies) made of stars and gas? 
Was each merger accompanied by strong star formation? Or 
conversely, do they originate from the early aggregation of 
lumps of gas turned into stars in the remote past via a burst- 
like episode ever since followed by quiescence so as to mimic 
a sort of monolithic process? If so, what is the key physical 
parameter governing this happening? Why do faint dwarf 
spheroidal and elliptical galaxies show evidence of recent 
star formation, whereas the luminous elliptical galaxies show 



The answers to the above questions may perhaps come 
from understanding the physical parameters governing the 
star formation history (SFH) of Dwarf Galaxies (DGs) of 
the Local Group. 



These objects in fact are thought to be the most com- 
mon type of galaxies in the Universe. Therefore, their stellar 
populations, their gas and dark matter (DM) content, their 
dynamical properties etc. are the subject of a great deal of 
observational and theoretical studies (see Mateo 1998 for a 
recent review on the subject) to unravel the mechanism of 
formation. 



From the theoretical side, much effort has been spent to 
infer the star formation history of DGs from their observed 
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CM-diagrams (CMD), and to this purpose different sophis- 
ticated techniques have been developed (Gallart et al. 1999, 
Hernandez et al. 2000 and references therein). 

Re-constructing the star formation history of these 
galaxies is of paramount importance in order to understand 
whether they formed in isolation or suffered from mergers or 
tidal stripping during their evolution. This galaxy types are 
indeed believed to be the building blocks for the assembly 
of larger galaxies. 

Recent investigations proved that in many cases the star 
formation history deduced from the CMD is very irregular. 
Nice examples are the dwarf spheroidal (DSph) galaxies Ca- 
rina and Leo I which together with a very old population 
show the signature of recent star forming activity. 

Other galaxies, however, seem to have experienced a 
single burst of star formation, although rather extended in 
time, like Draco, Ursa Minor or Sculptor, whose CMDs re- 
semble those of globular clusters in the Milky Way halo. 

It is often argued that DSGs were born almost at the 
same time, but evolved along different paths, probably de- 
pending on their distance to the Milky Way, their orbit and 
the amount of interaction experienced in the past (see Mateo 
1998). Is this the sole possibility? 

A sample of DSph galaxies is presented in Table [j], to- 
gether with a compilation of their basic properties. We show 
the CMD of two of them, Carina and Sculptor, in Fig. [jj 
They represent two very different evolutionary paths: in the 
case of Carina, star formation seems to have lasted more 
than about 6 Gyr, where in the case of Sculptor a dominant 
old burst of short duration seems to be more appropriate. 

In this paper, instead of analyzing CMDs, we propose a 
completely different approach, i.e. we intend to derive from 
first principles the star formation history of dwarf galaxies 
together with its large range of types. Starting from "real- 
istic" conditions for cosmological DM halos, we follow by 
means of N-body Tree-SPH simulations the formation of 
dwarf galaxies out of the cooling gas in the DM halo po- 
tential well, and the transformation of gas into stars. 

The present analysis stems from the recent study by 
Carraro et al. (2001), in which monolithic models of ellip- 
tical galaxies of different mass have been investigated by 
means of similar N-body Tree-SPH models with complete 
description of star formation, chemical enrichment, and en- 
ergy feed-back. The key result of those simulations is that 
starting from virialized haloes of DM with the same initial 
density, at decreasing galaxy mass the star formation his- 
tory gradually changes from a single sharply peaked initial 
episode to a series of burst stretching over most of the Hub- 
ble time. This is a very promising result because it consid- 
erably alleviates several points of weakness of SDGW model 
of Larson (1974) currently in use to interpret the data of 
elliptical galaxies. 

Applying the Carraro et al. (2001) scheme to the family 
of dwarf galaxies in which a large variety of star formation 
histories are seen for objects spanning a rather small range 
of masses, another dimension should be added to the prob- 
lem. This turns out to be the initial density. Combining the 
results of Carraro et al. (2001) with those we are going to 
present here, the monolithic collapse of baryonic material 
may lead to different star formation histories ranging from 
an initial prominent spike, to a sort of broad peak, to a se- 
ries of bursts of different intensity, and finally to an ever 



Figure 1. The Color-Magnitude diagram of two DSGs: on the left 
Carina, on the right Sculptor. This figure is taken from Smcckcr- 
Hane & McWilliam (1999). 



Table 1. Properties of some Dwarf Spheroidal Galaxies in the 
Local Group. Data from Mateo (1998). 



Galaxy 


M* 


r R 


<7 


[Fe/H] 




1O 6 M 


arcmin 


km/sec 




Carina 


13 


8.8±1.2 


6.8±1.6 


-2.0±0.2 


Phoenix 


33 






-1.9±0.1 


Leo I 


22 


3.3±0.3 


8.8±0.9 


-1.5±0.4 


Sextans 


19 


3.2 


6.6±0.7 


-1.7±0.2 


Ursa Minor 


23 


15.8±1.2 


9.3±1.8 


-2.2±0.1 


Draco 


22 


9.0±0.7 


9.5±1.6 


-2.0±0.15 


Sculptor 


7 


5.8±1.6 


6.6±0.7 


-1.8±0.1 



continuing, fluctuating mode. The above trend holds for a 
given initial density at decreasing galaxy mass, or for a given 
galaxy mass at decreasing initial density. 

The layout of the paper is as follows. Section 2 describes 
the numerical tools and the set-up of the initial conditions of 
the galaxy models. Section 3 presents the evolution and final 
properties of the galaxy models. Section 4 deals with the 
analysis of the CMDs that based on the star formation and 
chemical enrichment histories are expected to be observed 
at the present time and compares them with those of a few 
proto-type objects. Finally Section 5 summarizes the results 
and draws some general conclusions. 



2 N-BODY, TREE-SPH MODELS 

The simulations presented here have been performed using 
the Tree-SPH code developed by Carraro at al. (1998) and 
Buonomo et al. (2000). In this code, the properties of the 
gaseous component are described by means of the Smoothed 
Particle Hydrodynamics (SPH) technique (Lucy 1977 and 
Gingold & Monaghan 1977), whereas the gravitational forces 
are computed by means of the hierarchical tree algorithm of 
Barnes & Hut (1986), using the tolerance parameter 9 — 0.8, 
expanding the tree nodes to quadrupole order, and adopting 
a Plummer softening parameter. 

In the SPH method, each particle represents a fluid el- 
ement whose position, velocity, energy, density etc. are fol- 
lowed in time and space. The properties of the fluid are 
locally estimated by an interpolation which involves the 
smoothing length h. In our code each particle possesses its 
own time and space varying smoothing length h, and evolves 
with its own time-step. This renders the code highly adap- 
tive and flexible, and suited to speed-up numerical calcu- 
lations. No other details on the N-body Tree-SPH code in 
usage are given here for the sake of brevity. They can be 
found in Carraro et al (1998) and Buonomo et al (2000). 

2.1 Physical input 

Star Formation. The star formation rate (SFR) follows the 
law 
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1 ' dt dt t a 



(1) 



Table 2. Initial conditions for the galaxy models. 



where c* and t g are the so-called dimensionless efficiency of 
star formation and the characteristic time scale, respectively. 
In these models we adopt c* = 1. This characteristic time- 
scale t g is chosen to be the maximum between the cooling, 
t coo i, and the free-fall, tff, time scales 



tcooi 



E 



\E H -K c /p g \ V4ttGp 

where E is the current thermal energy per unit mass of gas, 
Eh is the heating rate per unit mass of gas by all possi- 
ble sources (supernova explosions, stellar winds, mechanical 
agents etc.), and Ac is the cooling rate (ergs -1 cm -3 ) by 
radiative processes. For more details on Eh, Ac and asso- 
ciated cooling time scale see Chiosi et al. (1998). Finally, 
the free-fall time scale is calculated by means of the total 
density p (gas, stars, and DM). Furthermore, a gas particle 
is eligible for star formation if the following conditions are 
met (Katz 1992; Navarro & White 1993): 



f— V 

\4ttG P J 



(2) 



tsound ^ tff 



and 



tcooi « tff 



where t BOUn d is the propagation time-scale of sound waves. 
It is worth recalling that t coo i, tff and t soun d ultimately 
depend (although in a different fashion) on the density of 
the fluid. 

Radiative cooling. The cooling rate is taken from the 
tabulations of Sutherland & Dopita (1994) and Hollenbach 
& McKee (1979) and amalgamated as in Chiosi et al. (1998). 
This allows us to account for the effects of variations in the 
metallicity of the fluid elements as a function of time and 
position. 

Energy Feed-Back. This originates from the SNae ex- 
plosions of type la and II (Greggio & Renzini 1983) and 
stellar winds from massive stars (Chiosi & Maeder 1986). 
To clarify how the energy feed-back is implemented in the 
code, we briefly discuss how a gas particle is changed into a 
star-particle. At any time step a gas particle eligible for star 
formation decreases its mass by the quantity 

c* At, 



Am 



1 — exp(- 



t„ 



(3) 



which simply follows from integrating equation (||) over the 
time step, and a new star-particle is created with mass 
m s = Am. This in turn is conceived as being actually 
made of a large number of smaller sub-units (the real stars) 
lumped together and distributed in mass according to a 
given initial mass function. In other words a star-particle 
is a single stellar population whose temporal evolution is 
well known (see for instance Bressan et al., 1994). For the 
purposes of the present study we adopt the initial mass func- 
tion of Salpeter (1955) over the mass interval from 0.1 to 120 
Mq. As time elapses, the real stars in any star-particle evolve 
and die, injecting gas (partially with the original composi- 
tion and partially enriched in metals) and energy by stellar 
winds and supernova explosions (both of type II and type I). 
Since the IMF is known these quantities can be easily cal- 
culated. Sharing of gas, metals and energy with neighbors is 
made according to the SPH formalism. While the new star- 
particle will now feel only the gravitational field, the left over 
gas-particle and all the others of the same type will continue 



Model 


M D 


M B 


Rd 


po 




1O 9 M 


10 9 M Q 


kpc 


Mq /kpc 3 


A 


0.9 


0.1 


4 


3.9 x 10 6 


Bl 


0.9 


0.1 


16 


5.8 x 10 4 


132 


0.9 


0.1 


26 


1.3 x 10 4 


B3 


0.9 


0.1 


35 


5.5 x 10 3 



to feel both gravitational and hydro-dynamical forces. The 
process can be repeated as long as there are gas-particles 
around eligible for star formation. 

Two comments are worthwhile here: (i) of the energy 
released to the ISM by a SN event (about 10 51 ergs), only a 
tenth of it is given to the inter-stellar medium and the rest 
is radiated away (see Thornton et al. 1998 for an exhaustive 
discussion of this topic), (ii) the energy injection by stellar 
winds may parallel that from SNae explosions. For all other 
details the reader should refer to the study by Chiosi et 
al (1998), whose prescriptions are adopted in the present 
models. 

Chemical enrichment. Supernova explosions and stellar 
winds increase the metallicity of the gas. The enrichment 
process is described by means of the closed-box model (Tins- 
ley 1980; Portinari et al. 1998) applied to each gas-particle 
of mass m g . In brief, each gas-particle is viewed as a close 
box system, in which star formation occurs turning gas into 
stars and increasing the metallicity according to the simple 
relation 

Z = y z ln[ mg ] (4) 
m g ~ m s 

where m v is the initial mass of particle (gas), m s is the 
mass of stars borne in it, yz is the so-called yield per stellar 
generation (constant as long as the IMF is constant). We 
adopt here y z = 0.004. See Tinsley (1980) and Portinari 
et al. (1998) for definitions, more details and recent esti- 
mates of yz ■ Sharing of the metals among the gas-particles 
is described by means of a diffusive scheme (Groom 1997; 
Carraro et al. 1998) . For more details on the above physical 
ingredients and their implementation in the numerical code 
see Buonomo et al. (2000 and references therein). 

2.2 Initial conditions 

The initial conditions of a galaxy are the total mass, Mt, 
made up in turn by DM and baryonic matter (BM) with 
masses Md and Mb respectively, the number of DM and 
BM particles, their positions and velocities, and the mean 
density of the system or equivalently its initial radius. All 
the simulations are made using 10,000 particles for both DM 
and BM. 

The following considerations are made to set up the 
initial conditions listed in Table ^. 

(i) We start from the moment in which a lump of DM 
and BM (in standard cosmological proportions) detaches it- 
self from the Hubble flow and begins to collapse toward the 
virial conditions. It is commonly assumed (see for instance 
Lacey & Cole 1993, 1994 and references therein) that this 
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occurs when the local density p[z) equals or exceeds the 
mean density of the Universe pu(z) by a certain factor 

*>0) = (Po) > 200 x p u (z) 

where z is the red-shift, (po) is the initial density of the 
proto-galaxy and pu(z) is 



3 h 2 x 100 2 3 

P"(2) = (1 + z ) 

Pu(z) = 1.99 x 10~ 29 h 2 {l + zf 



(5) 



g cm 



where ft = Ho /WO is the normalized Hubble constant. All 
other symbols have their usual meaning. 

(ii) Upon assigning the total mass Mt of the proto- 
galaxy, we may derive the initial radius, otherwise known as 
the virial radius or -R200, 



fJLws v ( Mt 
v 4tt ; V 200p u * 



a/3 



x (l + z)~ 



(6) 



mi / Mt 

#200(2) = 0.09617 x 200" 1/3 X I -j± 

where Mt is in solar units, and R200 is in kpc. In the fol- 
lowing we will adopt Hq= 65 km/s/Mpc and Zf or = 5. 

(iii) In order to simplify things we assume that the ini- 
tial lump of DM very soon acquires the condition of virial 
equilibrium, so that the distribution of DM in the proto- 
galaxy has the density profile 



p(r) 



pc- 



(7) 



where r c and p c are the radius and density, respectively, of a 
small central sphere. The spatial coordinates of DM particles 
inside the initial sphere of radius -R200 are assigned by means 
of a Monte-Carlo procedure. 

The initial velocity of a DM particle located at any ra- 
dial distance r is derived from the velocity dispersion o~(r) for 
a spherical isotropic collision-less system with the adopted 
density profile (Binney & Tremaine, 1987) 



p(r)a(r) = 



GM(r') , , 
P5 P( r ) dr 



Inserting equation (Q) we get 



cr(r) = p c r c Gr In 



(8) 



(9) 



where Rt is the initial total radius, i.e. feoo- Finally, the 
velocity v(r) is set equal to v(r) = |<r(r) assuming equipar- 
tition among the three components v x , v y , and v z . 

(iv) The BM inside the DM halo follows the evolution 
of the latter, however, no star formation is allowed to occur 
before DM has reached the virial condition. The particles of 
BM (gas at the beginning) are distributed homogeneously 
inside the DM halo with zero velocity field, thus mimicking 
the infall of primordial gas into the potential well of DM 
(White & Rees 1978). 

(v) The initial values for the mass of BM and DM are 
in the fixed proportion Mb — 0.1 Mo- Since the numerical 
simulations are made using the same number of particles 
(10,000) for both DM and BM, the mass resolution of the 
DM is 10 times lower. The initial metallicity of the gas is 
assumed to be Z = 10 -4 , the lowest value at our disposal in 
the grid of stellar models and related chemical yields. 

(vi) The initial total mass of our models is chosen in 



such a way that the present day total mass in stars agrees 
with the observational estimates for many dwarf galaxies of 
the Local Group. According to Mateo (1998; see also the en- 
tries of Table [l]) the mass in stars is about 10 to 30 xlO 6 M0. 
Basing on this we adopt the typical values Mt = 1 x 1O 9 M0, 
M D = 0.9 x 1O 9 M and M B = 0.1 x 10 9 Af Q . Owing to 
exploratory nature of this study the above choice is fully 
adequate to our purposes. 

(vii) There is a final consideration to be made con- 
cerning the initial density (radii) of our models. In prin- 
ciple, there should be one-to-one correspondence between 
red-shift and over-density of the proto-galaxies with respect 
to the surrounding medium. However, fluctuations around 
the mean excess are always present. For the most massive 
galaxies the degree of freedom is small because the fluctua- 
tions are small at the large scales. This is not the case for 
galaxies of lower mass because the fluctuations are larger at 
the smaller scales. From a technical point of view we will 
consider two classes of models: (a) a first group (labelled 
A) whose initial density strictly follows relation (H); (b) a 
second group (labelled B), in which the initial density is ar- 
bitrarily set to the value p„(.z)/200 and let fluctuate around 
it. Their initial radii are accordingly expressed by relation 
([]) in which the factor 200~ 1,/3 is dropped. It is worth stress- 
ing that this second group of models is only meant to explore 
the consequences of varying the initial density of the proto- 
galaxy. The analysis below will show that the initial density 
bears very much on the past SFH passing from a single ini- 
tial episode for high densities to recurrent burst-like episodes 
for low initial densities. 



3 THE MODEL GALAXIES 

The relevant data characterizing the final state of evolution 
of the four galaxies are summarized in Table ^. These are 
the mass in stars, the mass of remaining gas, the present- 
day effective radius of the baryonic matter R b ,b, the cen- 
tral velocity dispersions of baryonic and dark matter, at 
and o d respectively, and an estimate of the maximum and 
mean metallicity weighted on the star formation rate ^(t), 
\og(Z/Zo)Max and {\og{Z / Zq))^ , respectively. In the dis- 
cussion below particular attention is given to models of type 
B because they present the largest variety of SFHs. The re- 
sults for the case A galaxies are very similar to those of case 
Bl. 



3.1 Star Formation and Metal Enrichment 

Inside the virialized halo of DM, gas cools and sinks toward 
the center of the gravitational potential well, and sooner or 
later forms stars according to the adopted prescription. 

The four panels of Fig. ^ show the SFR of our galaxies 
as a function of time. In model A with the highest initial 
density the SFR reduces to a prominent, initial burst of 
activity followed by quiescence. All the action is within a 
narrow age range of about 0.2 Gyr. 

Passing to type B models the situation is more complex. 
Model Bl despite its significantly lower initial density com- 
pared to case A, still resembles this latter because the star 
forming activity is confined in the past. However, two im- 
portant new features appear: first star formation is delayed 
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Table 3. Properties of the Galaxy models. 



Model 


M a tar 


Mgas 






0"D 


Pc,B 


log(Z/Z B ) Max 


log(Z/Z Q )y 




10 7 M Q 


10 7 Af Q 


Kpc 


km/s 


km/s 


M /pc 3 






A 


3.2 


6.8 


0.07 


10 


25 




-0.017 


-1.158 


Bl 


7.7 


2.3 


0.24 


4.7 


18.7 


0.70 


-0.128 


-0.577 


B2 


2.5 


7.5 


1.55 


6.2 


15.2 


0.09 


-0.671 


-1.079 


B3 


1.3 


8.7 


8.70 


3.1 


12.7 


0.01 


-1.032 


-1.408 




4 2 
Look-Back Time (Gyr) 



Figure 2. The SFR history in the four galaxy models with dif- 
ferent mean initial density as a function of the look-back time in 
Gyr. The SFR is given in M0/year. The initial density decreases 
from the bottom to the top panels (cf . Table ^| ) . 



for about 1 Gyr and second it stretches over an age range of 
about 3 Gyr. The SFR is somewhat irregular and intermit- 
tent. Moving to even lower initial densities (models B2 and 
B3) the features in question are amplified, the SFR shows 
indeed many burst-like episodes of short duration (shorter 
than 0.1-0.2 Gyr) superimposed on a gentle, nearly periodic 
trend with typical time scale of about 1 Gyr. This behavior 
is natural in small size galaxies (Hirashita et al. 2000) and is 
the consequence of the balance between cooling and heating 
as a function of the depth of the potential well. More precise- 
ley, the rate of star formation and heating in turn scale with 




4 2 
Look-Back Time (Gyr) 



Figure 3. Metal enrichment of the three model galaxies of type 
B as a function of the look-back time. The filled square along each 
curve shows the mean metallicity weighted on the star formation 
rate. 



the density, whereas the cooling rate scales with the square 
of the density. Therefore, in a high density environment not 
only more stars are made and more energy and metals are 
injected but also cooling is more efficient so that it easier 
for the gas to meet the star formation condition. Star for- 
mation is somehow forced to completion. In contrast, in a 
low density medium, the rate of star formation and energy 
and metals injection tend to be lower, but also cooling be- 
comes less efficient. The gas particles find it more difficult to 
meet the conditions required for star formation. As a con- 
sequence of it, the star forming activity become fluctuating 
and even discontinuous. 

In the course of evolution, the galaxies convert their ini- 
tial gas content into stars at different efficiencies. In general, 
the higher the mean initial density, the larger is the fraction 
of gas converted into stars. However, this trend can be coun- 
teracted by the effects of heating by supernova explosions, 
stellar winds etc. which in the case of very high efficiency 
such as in our case A may even lead to a lower fraction of 
gas being converted into stars (see the entries of Table |^ to 
be described below). 
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-27 - 



13 3 4 

R (kpc) 



Figure 4. The final distribution of the star density (in g cm -3 ) 
in model B2. Superimposed is a Hernquist profile for the effective 
radius R e g = 1.55 kpc. 



In all cases under consideration, the gas left over by the 
star forming process is partly lost (out of the virial radius) 
and partly locked in the outer DM halo. Finally the left-over 
gas is partly enriched (30%) and partly unprocessed (70%) 
because of the partial diffusive mixing. 

In Fig. |^ we show the history of metal enrichment for 
the models Bl, B2 and B3. Model A is not displayed here 
because thanks to the very initial spike of stellar activity and 
quiescence ever since, the metallicity is practically constant 
with time. In all other models we notice a continuous but 
irregular increase, characterized by some pauses and sev- 
eral jumps, which mark successive bursts of star formation. 
The maximum metallicity goes from \og(Z / Zq) Max=-1 -032 
(B3) to log(Z/Z e ) Max =-0.128 (Bl) and log(Z/Z Q ) Max =- 
0.017 for model A. However, the mean metallicity weighted 
on the star formation rate (log(Z/ZQ))xn is much lower. It 
varies from (log(Z/Z Q ))*=-1.48 (B3) to (log(Z/Z©)>*=- 
0.577 (Bl) and (log(Z/Z©)}*=-1.158 (A). 

It is worth commenting here that these values of the 
metallicity are merely indicative of the real situation owing 
to our over-simplified treatment of chemical enrichment. 



3.2 Final structure of the model galaxies 

To illustrate how the mass of three main components, i.e. 
DM, gas and stars, is distributed in our model galaxies as a 
function of the radial distance, in Fig. |we show the present- 
day radial profiles of the mass density in gas (solid line), 
stars (dotted line), and DM (dashed line) limited to the 
cases of models Bl and B3. 
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Figure 5. Final radial distribution of the gas (solid line), stars 
(dotted line) and DM (dashed line) densities (in g cm -3 ) in the 
model galaxies Bl and B3. 

3.2.1 Dark Matter 

In our simulations, the DM halo surrounding the stellar and 
gaseous matter does not significantly change its properties 
in the course of evolution. The density profile of DM is very 
smooth, and depending on the initial conditions its central 
value ranges from a few 10 -25 g/cm 3 to a few 10" 23 g/cm 3 . 
The most important result is that the innermost regions of 
the galaxies (within about 1 kpc of the center ) are dom- 
inated by baryons, whereas the outer parts are dominated 
by DM (see Fig. ^) . Moreover, the velocity dispersion of the 
DM is about twice the velocity dispersion of the stars. These 
results agree with the notion that present-day DGs are DM- 
dominated (Salucci & Persic 1998), because the original con- 
tent of baryonic material has been largely lost in galactic 
winds (see below). 

3.2.2 Stars 

The final stellar mass of the models is 3.2 x 10 7 M for case A 
and it ranges from 1.3 x 10 7 Mq to about 5.5 x 1O 7 M0 going 
from Bl to B3. The stellar component always dominates in 
the innermost regions of the galaxy models. In order to check 
whether our models may resemble real dwarf galaxies, the 
dwarf spheroidals in particular, we look closely at the model 
B2. At the end of the simulation the central density of stars 
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Figure 6. Iso-density contours of the gas distribution in model 
B2 as seen at the end of the simulation. The X and Y-axis are the 
spatial coordinates in kpc centered on the center of the gravita- 
tional potential. Note the circular annulus void of gas surrounding 
the central body of the galaxy. The gas in the outermost region is 
beyond the virial radius and therefore escaping from the galaxy 
(galactic wind). 

is about 3 x 10~ 24 g/cm 3 and the stellar matter occupies a 
central sphere of about 5 kpc radius (see Fig. ||), which is 
about five times smaller than the initial virial radius, even 
though the density gradient is so steep that most of the stel- 
lar mass is within the effective radius R e ,B- The spherically 
averaged stellar density profile of the same model is then 
compared with the Hernquist (1990) profile for the stellar 
distribution in a dwarf spheroidal with effective radius 7? e ,s 
= 1.55 kpc (see Fig. Q). The remarkable agreement between 
the two profiles shows that the final galaxy resembles in- 
deed a dwarf spheroidal. The effective radius of the model 
is about 17 times less than its initial virial radius. Finally, 
we like to note that as for kinematics the stellar velocity 
dispersion in the center of our model galaxies ranges from 
10 to 3 km/s passing from case A to cases Bl, B2 and B3 
(see Table |). 

3.2.3 Gas 

Most of the initial gas originally inside the DM halo is not 
processed into stars (see the entries of Table . 

As expected in all models but for minor differences the 
spatial distribution of the gas in the central region closely 
follows the distribution of the stellar content (see Fig. ^) . 

However, there is an interesting feature to note looking 
at the large scale distribution of gas. The case is illustrated 
in Fig. |^ for the proto-type model B2. Of the gas left over 
inside the DM halo, only the fraction located in the central 
region keeps the original spatial distribution. The remaining 
part, heated up by the energy feed-back from supernova ex- 



Figure 7. Present-day gas and velocity profiles in the model 
galaxy B2. Upper Panel: The fractionary mass distribution of gas 
as function of the radial distance in kpc. Bottom panel: the radial 
velocity component V r (r) (dots) and the escape velocity V eS c(r) 
(long dashed line) as a function of the galacto-centric distance. 
The vertical dashed line is the virial radius, whereas the vertical 
solid lines show the effective radius of baryonic matter, i? e ,s and 
the layer, r esc , at which V r (r) > V eS c- All the gas located at 
T > Tesc freely flows out. 



Table 4. Galactic winds. 



Model 


Mg 


Re,B 


R200 


, e s c 


AMg tW 


AM,,. 
M g 




10 7 Af Q 


kpc 


kpc 


10 7 M o 


1O 7 M 


% 


A 


6.8 


0.07 


4 


2.4 


4.5 


66 


Bl 
B2 
B3 


2.3 
7.5 
8.7 


0.24 
1.55 
8.70 


16 
26 
35 


1.5 
5.2 
6.4 


0.8 
2.3 
2.1 


35 
30 
21 



plosions and stellar winds, reaches regions very far from the 
central sphere containing most of the stars (i.e. the luminous 
part of the galaxy). 

Finally, it is worth noticing that between the central 
part and the most distant regions there is a very low den- 
sity, almost spherical shell, which is due to the strong galac- 
tic wind episode occurred in the past. This gas distribution 
resembles indeed the distribution of cold gas recently dis- 
covered by Blitz & Robishaw (2000) in a sample of dwarf 
spheroidal galaxies of the Local Group. 

The major question is to be addressed whether or not 
part of the left over gas becomes gravitationally unbound to 
the main body of the galaxy thus freely escaping as galactic 
wind. 
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Table 5. T G is the present age (in Gyr) of the Galaxy models. 
Tjj,zfor is the age of the Universe at the epoch of galaxy for- 
mation. Hq is the Hubble constant in km/s/Mpc, whereas 50 is 
the deceleration parameter. Finally Zf or is the red-shift at which 
galaxies are assumed to form. 



H 


90 


z f or 


T G 


T U,z for 


50 


0. 


5 


16.450 


3.290 


50 


0.5 


5 


12.265 


0.895 


60 


0. 


5 


13.708 


2.742 


60 


0.5 


5 


10.220 


0.746 


70 


0. 


5 


11.750 


2.350 


70 


0.5 


5 


8.760 


0.640 



Figure 8. A simulated CMD according to the SFR of case A pre- 
sented in Fig. ^| The mean metallicity is (log(Z/Z0)}$=-1.158. 
The left panel is for H =70, <jo=0.5, z for =5, T G =8.76 Gyr and 
Tu,zior =0-64 Gyr. The right panel is the same but for Ho=50, 
q =0,°z/ or =5, T G =16.45 Gyr and T UiZ r = 3.29 Gyr. Finally, la 
photometric errors are displayed at the right side of each panel. 



3.2.4 Galactic winds 

The answer to the above question is given by Fig. [?], which 
displays the case of model B2, and Table ^, which lists the 
relevant data for all the models. 

The top panel of Fig. |?] shows the fractionary gas con- 
tent as a function of the radial distance in kpc, whereas 
the bottom panel shows the radial component of the gas 
velocity V r (r) (dots) and the escape velocity V eS c{r) (long 
dashed lines) as a function of the galacto-centric distance. 
The vertical, dashed line is the initial radius R200, whereas 
the vertical solid lines labelled R e B and V r > V esc show the 
effective radius of BM and the layer at which V r > V esc . 

There is a layer R e3 c at which the condition V r > V e3C 
is met. All the gas lying above it can freely leave the galaxy 
in form of galactic wind. However, the amounts of gas lost 
in the wind varies with the type of model i.e. the initial 
density. The situation is summarized in Table ^ which lists 
the model identification label, the total among of left-over 
gas M g , the effective radius of BM i? e ,s, the initial virial 
radius -R200, the gas content M g ^ ac up to at the layer at 
which the velocity equals the escape velocity, the amount 
of gas lost in the wind AMg lW , and the percentage of gas 
lost in the wind with the respect to the total available gas, 

In model A, the gas mass blown away in the galactic 
wind amounts to about 66% of the remaining gas, whereas 
in models of type B it is roughly a factor of two less. 

There is a final remark to be made: within the effec- 
tive radius of BM, very little gas remains. This is nearly 
zero for model A, it amounts to ~ 5% for model Bl, to 
~ 15% for model B2, and ~ 20% for model B3. If our 
case A can be taken as representative of dwarf spheroidal 
galaxies, this would be consistent with the observed lack 
of gas in these systems (Mateo 1998). In contrast the rel- 
tively higher gas content in models B could make them 
more suited to interpret dwarf irregulars (case B3) or dwarf 
spheroidals/irregulars (cases Bl and/or B2) that seem to be 
gas-rich (Mateo 1998), 



Figure 9. The same as ing, but for case Bl. The mean metallicity 
is (log(Z/Z ))tf=-O.577. The left panel is for iT =70, <jo=0.5, 
z for =5, T G =8.76 Gyr and T UtZfor =0.64 Gyr. The right panel 
is the same but for Ho=50, qo=0, Zf or =5, T G =16.45 Gyr and 
T Ufor = 3.29 Gyr. 



Figure 10. The same as ing, but for case B2. The mean metallic- 
ity is (log(Z/Z ))*=-l.O79. The left panel is for H o =70, <?o=0.5, 
z for =5, T G =8.76 Gyr and T v>z r =0.64 Gyr. The right panel 
is the same but for Ho=50, <?o = 0, Zf or =5, T G =16.45 Gyr and 
Tu,z for = 3.29 Gyr. 

4 COLOR-MAGNITUDE DIAGRAMS 

The aim of this section is to check whether the SFR and en- 
richment law of our model galaxies may generate CMDs like 
those observed today in dwarf galaxies of the Local Group. 

Prior to this, one has to choose the link between the 
red-shift of galaxy formation, zf or and age in order to con- 
vert the rest-frame ages used so far into absolute ages of 
the stars which would now be observed. This means that a 
certain model for the Universe has to be adopted. For the 
purposes of the present study we adopt the Friedman model 
described by the Hubble constant Ho and deceleration pa- 
rameter qo- In all simulations we assume the same red-shift 
of galaxy formation, i.e. zj or — 5. To facilitate the use of 
the simulations below, in Table |^ we summarize the maxi- 
mum ages T G of the galaxy models, and the age Tu,z for of 
the Universe at Zf or = 5, for three different values of Ho 
(50, 60 and 70 km/s/Mpc), and two values of qo (0 and 0.5). 
Higher values of both quantities are not likely here because 
they would predict galactic ages T G that are too old, at least 
compared to what is currently known about the age of glob- 
ular clusters, i.e. from 10.0 to 15.8 Gyr as recently reviewed 
by Carretta et al. (2000). 

Synthetic CMDs are constructed by means of a popula- 
tion synthesis algorithm. Given the galaxy age To, for each 
value of the stellar age t, an isochrone with the mean metal- 
licity (Z) is derived from a set of stellar evolutionary tracks. 
These are from Girardi et al. (2000) for 0.0004 < Z < 0.03, 
and Girardi et al. (1996) for Z = 0.0001. The evolutionary 
sequences span the mass and age ranges of 0.15 Mq < M < 
7 Mq, and from 0.06 to 20 Gyr. The mass and age ranges are 
fully adequate to describe the stellar populations of dwarf 
spheroidal galaxies. All important evolutionary phases are 
included, from the zero-age main sequence to the end of the 
thermally-pulsing asymptotic giant branch. The isochrones 
in the theoretical plane (luminosity versus effective temper- 
ature) are transformed to the observational magnitudes and 
colors using the transformations provided by Bertelli et al. 
(1994). 

Each isochrone of age t is populated by stars according 
to the Salpeter initial mass function (IMF), and in propor- 
tions fixed by the star formation rate, SFR(T G — t). The rel- 

Figure 11. The same as in^, but for case B3. The mean metallic- 
ity is (log(Z/Z ))*=-1.4O8. The left panel is for H =70, <?o=0.5, 
z for =5, T G =8.76 Gyr and is T [/jZ/or ,=0.64 Gyr. The right panel 
is the same but for Ho=50, qo=0, Zf or =5, T G =16.45 Gyr and 
Tu,z for = 3.29 Gyr. 
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ative number of stars predicted in each section of the CMD 
is added to the composite population. The procedure con- 
tinues for all values of the age, i.e. from the present (t = 0) 
up to the maximum galaxy age (t = To). 

Finally, the CMD is populated using a arbitrary total 
number of stars. In the present simulations, we include the 
effects of photometric errors and neglect those of binaries for 
the sake of simplicity. Photometric errors are assumed to be 
similar to those obtained by Smecker-Hane & McWilliam 
(1999) in HST observations of Carina dSph. 

The most uncertain ingredient of these CMD simula- 
tions is mass loss along the red giant branch (RGB) of low- 
mass stars. This is taken into account by decreasing the 
stellar mass when passing from the tip of the RGB to the 
zero-age horizontal branch (see Renzini 1977). Mass loss is 
assumed to follow the Reimers (1975) law with a suitable 
value of the parameter r\ (Fusi-Pecci & Renzini 1976). It is 
worth recalling that by applying mass loss at the stage of 
isochrone construction (and not during the computation of 
evolutionary tracks) allows us to easily vary r\ from simu- 
lation to simulation if required. We adopt the mean value 
(77) = 0.50 with dispersion a v = 0.05. These values are 
slightly larger that those used to reproduce the morphol- 
ogy of globular clusters (Renzini & Fusi-Pecci (1988) but 
still within current uncertainties. 

The above technique is used to generate the CMDs 
shown in Figs. | through |ll] for the SFHs (Fig. |) and the 
enrichment rate (Fig H) of our model galaxies A, Bl, B2, 
and B3. Furthermore, CMDs are presented for two, some- 
what extreme values of the galaxy age Tg corresponding to 
different values of Ho, qo, Zfor, Tu,zfor as indicated. 

Model A. Case A is straightforward because owing to 
its prominent and unique initial burst of activity, the associ- 
ated CMD would simply resemble that of a globular cluster. 
This is shown in the two panels of Fig.^. The left panel is 
for _ffo=70, qo=0.5, Zf or =5 and Tg=8.76 Gyr, whereas the 
right panel is for Ho=50, qo=0, Zf or =5, and Tc=16.45 Gyr. 
The two panels clearly show how the CMD of a globular 
cluster with intermediate mean metallicity ((log^/Z0)*=- 
1.158) changes as a function of the age. Passing from the 
young (left) to the old (right) population, the turn-off lumi- 
nosity decreases whereas the clump of red stars (He-burners) 
turns into an extended horizontal branch. This is more pop- 
ulated at the red side because the metallicity is somewhat 
higher than the typical value in globular clusters with very 
extended horizontal branches, mostly populated at the blue 
side. 

Model Bl. The left and right panels of Fig.^| show the 
expected CMD for the models Bl. The CMD of the older 
age (right panel) much resembles the one of an old globular 
cluster with relatively high metallicity, (log Z/Zq)^ =-0.577 
according to the entries of Table |§| There is some blurring 
at the turn-off and sub-giant branch due to the large age 
spread caused by the prolonged initial stellar activity (about 
3 Gyr). These effect adds a new dimension to the horizontal 
branch, because it is now populated by stars with signif- 
icantly different initial mass, which tends to evenly popu- 
late the horizontal branch. Despite the high metallicity, this 
CMD could closely mimic the ones of Sculptor (see Fig. [j]) 
or Ursa Minor and Draco, i.e. galaxies dominated by a single 
burst of SF which extended for some time (see Hernandez 
et al. 2000). By decreasing the age down to 9 Gyr, the red 



clump (typical of an intermediate age population) appears. 
Although the star formation lasted for about 3 Gyr, its ever 
decreasing efficiency past the initial peak does not allow the 
presence of a secondary distinct generation of stars but only 
a large scatter at the turn-off. This CMD is a good candi- 
date to explain the properties of Leo II (Hernandez et al. 
2000). 

Model B2. The panels of Fig.|l(] show the expected CMD 
for the same value of zj or , Ho and qo=G as in Figs. |§|and ^ 
but for model B2. The CMD of the right panel does not sig- 
nificantly differ from its analog in Fig. ^ but for a larger scat- 
ter at the turn-off, sub-giant branch, and horizontal branch. 
The CMD in the left panel owing to its much younger age 
clearly shows the several distinct bursts of star formation 
already noticed in Fig. |^. 

Model B3. Finally we have the case of models B3, whose 
CMD is shown in the left and right panels of |ll] for the 
same combinations of Ho, qo and Zf or as in the previous 
simulations. The results of the right panel (_ffo=50, qo=0, 
Tq = 16.45 Gyr) are particularly interesting because three 
distinct old populations are now visible going from old to 
relatively young (say about 6 Gyr). Galaxies like Leo I or 
Carina could be reproduced by a star formation history like 
this. See for instance the CMD of Fig. [l] for Carina, or that 
of Leo I given by Held et al. (2000). The CMD of the left 
panel, owing to the much younger ages involved and almost 
continuous star formation up to the present could perhaps 
mimic the case of a DIrr. 

A general remark. The above simulations of dwarf 
galaxies and associated CMDs were not designed to match 
any particular dwarf galaxy. In most cases, the metallicity 
of the models was too high compared to observational value. 
This is certainly due to the adoption of the closed-box de- 
scription. Work is in progress to improve upon the treat- 
ment of chemical enrichment. Our simulations are intended 
only to show that dwarf galaxies with the same total mass of 
1O 9 M0 but originated from proto-galaxies of different initial 
densities may undergo different histories of star formation 
which in turn yield different CMDs, in agreement with the 
large variety of observational data for such galaxies. Dwarf 
galaxies formed in a high density environment tend to have 
a single initial burst of stellar activity. In contrast, galaxies 
of the same mass but with low initial density tend to un- 
dergo star formation over long periods of time, showing a 
number of burst-like episodes. In some extreme cases star 
formation may stretch up to the present. Accordingly we 
pass from galaxies that would nowadays appear as DSph, 
DE and DIrr. 



5 SUMMARY AND CONCLUSIONS 

We have tried to understand the physical reason for the 
large variety of star formation histories indicated by the 
stellar content of dwarf galaxies, with particular attention 
to the case of dwarf spheroidals and ellipticals. Indeed ob- 
jects whose estimated mass in stars, mean metallicity, and 
velocity dispersion span a rather narrow range of values, can 
have very different star formation histories. 

With the aid of Tree-SPH, N-Body simulations of 
galaxy formation and evolution we find that systems of the 
same total (dark and baryonic) mass but different initial 
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densities may have different star formation histories rang- 
ing from a single initial episode to prolonged stellar activity 
extended into a series of bursts over a large fraction of the 
Hubble time. 

Starting from "realistic" initial conditions, we follow the 
monolithic collapse of baryons inside the virialized halo of 
DM, and the formation of stars by gas cooling and collapse. 
The models include accurate treatments of star formation, 
chemical enrichment, energy feed-back from supernova ex- 
plosions and stellar winds, and metal-dependent cooling. 

The structural properties of the models remarkably 
agree with the observational data of real galaxies (mass den- 
sity profiles, central velocity dispersions, mean metallicities, 
mass-to light ratios, and even gas distributions). 

Furthermore, using the theoretical star formation and 
metal enrichment histories of the model galaxies we simu- 
late the CMDs that one would expect to observe today as a 
function of the underlying cosmological parameters Ho, qo, 
Zfar an d To - In many cases the simulated CMDs remarkably 
resemble the observational ones. 

The main conclusion of this study is that the adopted 
scheme of galaxy formation, i.e. early virialization of DM 
haloes and monolithic collapse of baryons into their potential 
wells, may lead to a variety of different situations as far 
as the detailed star formation histories is concerned, that 
depend on the initial density of the system. 

Together with the results of the companion paper by 
Carraro et al. (2001), briefly summarized in Section 2, this 
study suggests that total mass and initial density are the 
key parameters governing the final stellar content of ellip- 
tical galaxies. Ample possibilities exist for star formation 
histories, ranging from a single initial episode of short du- 
ration, to large time intervals, to burst-like activities, and 
even to very prolonged star formation. 

A final consideration is worth being made here. 
Throughout the analysis, when referring to the initial den- 
sity this has been linked to red-shift via the over-density with 
respect to the mean density of the Universe required by a 
proto-galaxy to detach himself from the Hubble flow and to 
collapse. This indeed was the over-density assigned to model 
A. Strictly speaking, at the red-shift of galaxy formation we 
have considered only models of type A are possible. To have 
models of type B brought into existence one has to wait for 
a fall-off of the over-density by a factor of 70 at least, or 
equivalently Azf or ~ 3. In other words, galaxies of similar 
mass but formed at zj or ~ 2 would behave as models B and 
for 50 < Ho < 60 and qo — would have a maximum age 
for their stellar populations still fully compatible with the 
information derived from the observational CMDs. In other 
words, proto-galaxies of the same total mass collapsing at 
different red-shifts would contain different kinds of stellar 
populations. 
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